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We have analyzed two long BeppoSAX observations of the bright Seyfert galaxy NGC 4151, searching for short 
timescale (~10-200 ksec) X-ray spectral variability. The light curve of a softness ratio, chosen as most sensitive to 
pinpoint changes of the column density of the absorbing gas along the line of sight, shows significant variations. 
We try to model these variations by performing a detailed, time resolved, spectral analysis. We find significant, 
large (factors of 1.5-6) variations of the absorber column densities on time scales of 40-200 ksec. These values are 
10-100 times shorter than those found by Risaliti et al. 2002 in a sample of Seyfert 2 galaxies, and provide strong 
constraints on the geometry of the obscuring medium. 



1. Introduction 

The geometry of the obscuring matter in Ac- 
tive Galactic Nuclei, predicted by unified models 
PP, is still unclear. The study of variability of 
the absorbing column density (N#) in these ob- 
jects, provides invaluable informations on the size 
and location of the absorber. Recently, Risaliti et 
al. [2] found Nh variability on time scales from 
~ 2 month to ~ 5 years in nearly all the objects 
belonging to a sample of bright, Compton-thin 
Seyfert 2 galaxies with multiple X-ray observa- 
tions. This variability suggests that the absorber 
could be much closer to the nucleus than usually 
assumed. 

To further investigate this issue, we have an- 
alyzed two BeppoSAX long observations of the 
brightest Seyfert galaxy, NGC 4151, and present 
here some preliminary results. 



2. Light curves analysis 

The two BeppoSAX observations were per- 
formed on 2001 December, 18 and 1996 July, 6 
with the Narrow Field Instruments: LECS (0.1- 
10 keV), MECS (1.3-10 keV), and PDS (13-200 
keV). 

The total absorbing column density mea- 
sured in NGC4151 by HEAOl, EXOSAT, Ginga, 
ROSAT and ASCA |4I5I6I7I8I91TU| ranges from a 
few 10 22 cm~ 2 to a few 10 23 cm~ 2 . Since the 
corresponding photoelectric cut-off lies in the 2-4 
keV energy range, we have chosen to analyze the 
light curves of the count rate in the bands 2-4 keV 
(MECS), 6-10 keV (MECS) and of their softness 
ratio (see Figures 1 and 2), to investigate whether 
Njj variations might have occurred within each 
observation. We also analyzed the light curve 
of the count rate in the 15-200 keV band (PDS) 
to monitor the continuum. From the analysis of 
these light curves we have selected some time in- 
tervals where we find indications that N# might 
have undergone significant changes. These inter- 
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vals are shown and labelled in figures 1 and 2. 
To quantify the spectral variations we have then 
performed a detailed spectral analysis of these 
stretches of data. 
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Figure 2. Same panels of figure 1 but for the 1996 
July 6 observation. 
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Figure 1. Light curves of 2001 December 18 
observation. From top to bottom: MECS 2-4 
keV count rate, MECS 6-10 keV count rate, 2- 
4 keV/6-10 keV softness ratio, PDS 15-200 keV 
count rate. The capital letters at the top of the 
figure, indicate the time intervals selected for a 
time-resolved spectral analysis. 



3. Spectral analysis 

Standard data reduction was performed using 
the SAXDAS software package version 2.0 fol- 
lowing Fiore, Guainazzi & Grandi [TT]. Spec- 
tral fits were performed using the XSPEC 11.2.0 
software package and public response matrices as 
from the 1999 December release. The PDS data 
were reduced using the "variable risetime thresh- 
old" technique to reject particle background [TT], 



LECS and MECS spectra were rebinned follow- 
ing two criteria: a) to sample the energy resolution 
of the detectors with three channels at all ener- 
gies whenever possible, and to obtain at least 20 
counts per energy channel. Constant factors have 
been introduced in the fitting models in order to 
take into account the intercalibration systematics 
between the instruments ^I] • The model adopted 
for the analysis has the following components: 

1. A power-law with an exponential high- 
energy cut-off; 

2. a neutral Compton reflection component 
(modelled with PEXRAV 

3. a narrow iron K a emission line with energy 
set to 6.38 keV in the observer frame; 

4. two components for the soft excess emis- 
sion: a thermal component, and a power 
law, reflection component. 
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5. two components for the intrinsic absorber, 
both assumed neutral: Nh ± which covers 
the nucleus totally; Njf 2 which covers the 
nucleus only in part, with a covering factor 
CV |4I5I6I7I8I9I10| . The absorber is com- 
pleted with the line of sight column den- 
sity through the Galaxy (N Hgal =2.17 X 10 20 
cm- 2 [Ej). 

Figures 3 and 4 illustrate the results for the 
two observations. In the 2001 December, 18 ob- 
servation, we find that from A to B to C, both 
absorbers and the photon index do not show sta- 
tistically significative changes. From C to M the 
two absorbers Njjj and N# 2 change at a confi- 
dence level ~ 90% by a factor of 2 and 3.3 re- 
spectively, while the photon index is constant at 
a confidence level ~ 99%. If the covering factor 
CV is kept constant throughout the fit the ampli- 
tude of the variations of Nh 1 and N# 2 is reduced 
to a factor 1.4 and 2 respectively. The variations 
occur on a timescale of 100-150 ksec. 

In the 1996 July, 6 observation, we find that 
from A to B both absorbers undergo changes at a 
confidence level ~ 90%, while from A to C the ab- 
sorbers change at a confidence level > 99%. From 
A to B to C, the photon index is again statisti- 
cally consistent with a constant. The amplitude 
of the variation of Nj^ and N# 2 is a factor 6 and 
5 respectively. If the covering factor CV is kept 
constant throughout the fit the amplitude of the 
variations of Ng^ and Njj 2 is reduced to a factor 
4 and 3 respectively. The variations occur on a 
timescale of ~ 30 ksec from A to B and sa 200 
ksec from A to C. 

4. Conclusion 

We have performed a detailed time-resolved 
spectral study of two long BeppoSAX observa- 
tions of the Seyfert galaxy NGC4151. We find 
strong evidences for spectral variability which we 
model in terms of variations of two absorbing 
screens. The timescales of these variations are be- 
tween 40 ksec and 200 ksec, 10-100 times shorter 
than the typical timescales found by Risaliti et 
al.[2] in a sample of Compton thin Seyfert 2 galax- 
ies. If we associate these timescales to a crossing 
time, the resulting linear size would be of about 
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Figure 3. From top to bottom: 68, 90, 99% \ 2 
confidence contours for N# 2 versus N/^ , N# 2 ver- 
sus photon index T and N/^ versus T, for the 
spectra corresponding to the intervals A, B, C 
and M in the 2001 December, 18 observation. 
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Figure 4. From top to bottom: 68, 90, 99% x 2 
confidence contours for N# 2 versus , Nh 2 ver- 
sus photon index T and N^ versus T, for the 
spectra corresponding to the intervals A, B and 
C in the 1996 July, 6, observation. 



one light day, which is the size of the innermost 
part of the Broad Line Region, as inferred from 
reverberation mapping in NGC 4151. This inter- 
pretation is supported by the argument, used by 
Risaliti et al. .2. , which assume that the absorb- 
ing medium might be made up of spherical clouds 
moving with Keplerian velocities around the cen- 
tral black hole. This requires that the cloud gas 
density is of the order of 10 9 cm~ 3 , typical of 
Broad Line Region clouds |13| . 
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